QCD at small baryon number 
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We consider the difficulties of finite density QCD from the canonical formalism. We present results for small 
baryon numbers, where the sign problem can be controlled, in particular by supplementing the pi=0 sampling with 
imaginary [i ensembles. We initiate the thermodynamic study of few-nucleon systems, starting with the measure- 
ment of the free energy of a few baryons in the confined and deconfined phase. We present a simple model for 
the phase transition, whose results are in good agreement with the literature, but extend to lower temperatures. 



1. INTRODUCTION 

QCD at finite density is of great physical inter- 
est, but numerically challenging: the usual Monte 
Carlo sampling weight becomes complex, causing 
a sign problem. Recently, progress has been made 
using different approaches |1I2I3I4| . However all 
these methods have severe limitations which are 
already being reached, with little prospect for fur- 
ther progress. 

Therefore, we propose yet another method - 
in a canonical formalism. Physically, we are in- 
terested in the study of systems consisting of a 
few baryons, for example the formation of nuclear 
matter or the binding of Deuteron. Here, we mea- 
sure the free energy of sectors with small baryon 
number, and present a simple model for the QCD 
deconfinement phase transition. Numerically, our 
sign problem is governed by the baryon number 
B, and not by the volume, and its dependence on 
temperature is mild. The sector B = (no sign 
problem) was analysed in detail in jo]. Here, as a 
continuation, we present results for B = 1, 2, . . . 

2. QCD'S PARTITION FUNCTIONS 

The grand canonical partition function with 
fixed chemical potential /i 

Z GC {p) = J[DU}er s ^> u] detM(U;n) (1) 

is numerically problematic for real chemical 
potential, since det M(U;n) becomes complex. 
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However, for imaginary /i = i/i/, det M(U; 
remains real. Moreover, Zac{i^i) is even and 
2^-periodic in fij 0. 

The canonical partition function with fixed 
quark number Q can be obtained from the grand 
canonical one using a Fourier transform, 

(2) 

The ^^-periodicity of Zcci^i) implies that 
Zc{Q 7^ mod 3) = 0. Hence, we only deal with 
integer baryon number B = Q/3. 

Knowing the canonical partition functions, 
ZccifJ-) can then be calculated for arbitrary 
chemical potential using the fugacity expansion 
(Laplace transformation): 

oo 

Z GC {n) = £ e- 3B *Z c (B). (3) 

B = -oo 

3. THE METHOD 

Our main observable is the free energy per 
baryon, F(B), defined by ^fra^o) = e~T F ( B ). 
We obtain Zc(B) by following the approach in 
jZj, which samples Zccii/J-i) at fixed fii's, then 
Fourier transforms each det(C7;i/ij) (all the ^in- 
dependency is in the determinant), 
v 

det M(U;ifj,i)= £ Z C (U; B)e- aB ^ . (4) 

B = -V 

Substituting into Eq. JTJ then into Eq. {2} gives 
Z C (B) = [[DU}Z c (U;B)e- s ^ . 
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We sample the grand-canonical ensemble at 
couplings [3 and ifii, hence the observables we 
measure are ratios of partition functions, 

Zc(P,B) = Z C (B,U) 

In the idea of the existence of a particu- 
lar fJ>i(B) to extract Zc{(3, B) optimally has been 
formulated. Here, in order to improve the accu- 
racy on the Zc(fl, By a, we simply combine results 
from several (/?, /Uj) -ensembles by Ferrenberg- 
Swendsen reweig hting|l. This requires knowing 
the determinant for different /i/'s. This informa- 
tion comes as a by-product of our determination 

of the Zc(B) in Eq. Q. In the temporal gauge 
(C/4(x, t) = 1 except for t — N t — 1), the staggered 
fermion matrix M can be written as 
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Following TO] , we define a fermionic transition 
matrix 

/ / ; i \ \ 

UNt-l 

with eigenvalues Ai,...,A6V- We obtain the 
Zc(U;B) using Eq. (0J and the equality 

det M(U; n) = e 3V % J[ (\i + e~$ 

4. RESULTS 1=1 

Since this is a feasibility study, we consider a 
small (4 4 ) lattice with volume V ~ (1.2 fm) 3 . 
We use four flavours of staggered quarks with 
mass m — 0.05a _1 . This theory is expected to 
have a first-order phase transition in the whole 
(T, n) plane. With Hybrid Monte Carlo, we gen- 
erate 55 ensembles of 4000-16000 configurations, 
at 11 temperatures y- G [0.8, . . . , 1.1], each at five 
imaginary chemical potentials. 

In the high-temperature phase T > T c , we ob- 
serve the characteristics of a system of free, rela- 
tively light quarks: the free energy F(fi) is pro- 
portional to /j, 2 for small /i J3J, so that the free 
energy F(B) ~ B 2 for small B (fi is conjugate 
to B). Furthermore, F(B) is almost indepen- 
dent of temperature in this regime. In Fig. 1, 
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Figure 1. Dimcnsionless free energy per baryon, 
T~ 1 F(B)/B, as a function of /? = 6/5 2 . 

we show T~ 1 F(B)/B as a function of (3 for vari- 
ous B, and indeed observe almost equidistant, p- 
independent levels (F{B)/B ~ _B). For entropy 
reasons, we expect = 1) to vanish in the 

thermodynamic limit. Deviations from these ex- 
pectations are caused by finite size effects. 

As T approaches T c from below, the free en- 
ergy drops dramatically. This is caused by an 
increase in the number of thermally accessible 
states, which triggers an entropic catastrophy a 
la Hagedorn. This collective phenomenon (of glu- 
ons) is visible even in the 1-baryon sector. 

In the low-temperature phase T < T c , our 
results are consistent with a simple model 
of weakly interacting heavy baryons at rest 
(Pmin = x a 1 ~ ja _1 ). At fixed temper- 
ature, F(B)/B stays constant over a whole range 
of B > 2, so that the free energy grows linearly 
with |B| as sketched in Fig. 2. The interaction is 
repulsive: F{B)/B - F(B = 1) » OTa" 1 - 60 
MeV per baryon, see Fig. 1. In the presence of 
a (real) quark chemical potential /i, the free en- 
ergy acquires an additional term 3/i-B, which tilts 
the free energy in Fig. 2. The free energy mini- 
mum jumps away from B = 0, and a first-order 
transition to a finite density occurs, when 

F(B). 



fir 



I plateau 



(5) 



The phase diagram resulting from our simple 
model is shown in Fig. 3. The agreement with 
other methods for /i < 0.35a _1 is very good, 
except at low /i, where our method suffers from 
finite size effects. In the thermodynamic limit, 
F(B = 1) will develop a singularity exactly at the 



3 



T>T„ 





T<T„ 





Figure 2. Sketch of the free energy F(B) in the 
two phases, for increasing (real) chemical poten- 
tial. For T < T c , the free energy minimum jumps 
to a non-zero density at a critical /i c . 



(T = T c , /-J. = 0) transition. Moreover, in the limit 
T -> 0, F(B = 1) = m N , yielding via Eq. © 
/i c = rnjv/3, which already is a good approxima- 
tion of the true phase transition. Therefore, we 
expect Eq. JSJ) to give a good estimate of the crit- 
ical chemical potential in the whole (T, /i) plane. 
On our small lattice, a plateau for F(B)/B is 
already reached for B = 2. The resulting phase 
boundary is shown in Fig. 3. Note the bending 
down for fi/T > 1, where the reliability of other 
methods may be questioned. 

5. CONCLUSION 

In 4-flavor QCD, we have observed that the 
free energy grows linearly with the baryon num- 
ber, for small densities at T < T c . This behav- 
ior is consistent with a first-order transition from 
weakly interacting heavy baryons to light, quasi- 
free quarks. The corresponding phase boundary 
agrees with the literature for moderate n/T, then 
bends downward. 

The same information could, in principle, 
be obtained by reconstructing the full grand- 
canonical partition function Eq.©. However, the 
sectors with large baryon number, which we did 
not consider, bring in little information but a lot 
of noise. This is why our approach works better. 

Since the sign problem is moderate, we can 
increase the volume and lower the temperature. 
As we do so, we expect the weak Pauli repulsion 
among baryons to be overcome by nuclear attrac- 
tion, leading to the formation of nuclear matter. 
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Figure 3. The (/x, T) phase diagram. An ap- 
proximation of the phase boundary is given by 
F(l)/3, where we neglect the baryonic interac- 
tions. F(2)/6 already gives the plateau- value in 
Eq. (J5J. Also shown are P and @|, which use 
the same simulation parameters. Agreement is 
very good for \i < 0.35a -1 . Ch. Schmidt re- 
cently obtained one point ^3] using the factoriza- 
tion method [T4*). 
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